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ABSTRACT 

In this paper, we present an expanding disc model to derive polarization properties of 
the Crab nebula. The distribution function of the plasma and the energy density of the 
magnetic field are prescribed as function of the distance from the pulsar by using the 
model by Kennel and Coroniti (1984) with a = 0.003, where a is the ratio of Poynting 
flux to the kinetic energy flux in the bulk motion just before the termination shock. 
Unlike previous models, we introduce disordered magnetic field, which is parameterized 
by the fractional energy density of the disordered component. Flow dynamics is not 
solved. The mean field is toroidal. 

Averaged polarization degree over the disc is obtained as a function of inclination 
angle and fractional energy density of the disordered magnetic field. It is found for 
the Crab that the disordered component has about 60 percent of the magnetic field 
energy. This value is also supported by the facts that the disc appears not 'lip-shape' 
but as 'rings' in the intensity map as was observed, and that the highest polarization 
degree of ~ 40 percent is reproduced for rings, being consistent with the observation. 

We suggest that because the disordered field contributes rather pressure than 
tension, the pinch force may be over-estimated in previous relativistic magnetohy- 
drodynamic simulations. Disruption of the toroidal magnetic field with alternating 
direction, which is proposed by Lyubarski (2003), may actually takes place. The rel- 
ativistic flow speed, which is indicated by the front-back contrast, can be detected in 
asymmetry in distributions of the position angle and depolarization. 

Key words: pulsars: general - ISM: individual: Crab nebula - supernova remnants 
- radiation mechanisms: non-thermal - polarization. 



1 INTRODUCTION 

It is well established that the Crab Nebula shines in syn- 
chrotron radiation. However the origin of the high energy 
particles has been a long-standing problem. The pulsar wind 
from the central pulsar is certainly the source of the parti- 
cles and the magnetic fields. We nevertheless do not know 
what fraction of the wind energy is in magnetic field and 
in plasma kinetic energy. We sometimes introduce the a- 
parameter, which is the ratio of the magnetic energy flux 
to the kinetic one of the wind bulk motion just before the 
termination shock. 

Kennel and Coroniti (KC) (1984a, 1984b) proposed a 
simple model in which the nebula is expanding hot plas- 
mas after the termination shock of the wind. In their model, 
the kinetic energy of the wind bulk motion converts into 
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heat with non-thermal components. This picture results in 
a small value of a, i.e., dominance of the kinetic energy, to 
explain the observed synchrotron luminosity and the pre- 
sumed expansion velocity of ^ 2000 km s~^. One finds it 
difficult to reproduce the kinetic energy dominated wind by 
means of the axisymmetric ideal-MHD model (Michel 1969, 
Bogovalov 1998, Begelman & Li 1994, Tomimatsu 1994), al- 
though some possibilities for the small a have been pointed 
out (Okamoto 2002, Vlahakis 2004). This is referred to as 
the cr-problem. 

Because the magnetic axis inclines to the rotation axis, 
the magnetic neutral sheet in the equator is folded for each 
rotation to produce a series of current sheets flowing out 
with the wind. Some authors attempt to solve the cr-problem 
by introducing dissipation of the magnetic field in such cur- 
rent sheets (Coroniti 1990, Kirk & Skjaeraasen 2003). 

Lyubarski (2003) has proposed that conversion of the 
magnetic energy to heat takes place when the current sheets 
come into the termination shock. With this picture, Poynt- 
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ing energy directly changes to plasma heat and non-thermal 
particles, so that the pulsar wind should not always be ki- 
netic energy dominant, i.e., the cr-problem may not be a 
problem. 

Magnetic field of the wind is essentially toroidal far be- 
yond the light cylinder owing to rotation. The nebula field 
convected from the wind is also believed to be toroidal, indi- 
cated by highly axisymmetric structures in the Crab Nebula. 
However, if disruption of the numerous current sheets takes 
place at or after the termination shock, one may expect that 
the nebula field is not pure toroidal but it has strong disor- 
dered components. Thus, the disordered magnetic field may 
be a good indicator for energy conversion from the toroidal 
magnetic field to plasmas. To diagnose disruptive processes 
of the magnetic field, it must be very much helpful to mea- 
sure magnitude of the disordered field. 

Disordered magnetic field can be detected by image 
analysis and by polarization analysis. Shibata et al. (2003; 
Paper I) calculated images by use of the KC model and com- 
pared with the Chandra observations. If the nebula field is 
pure toroidal, intensity along the major axis of the torus 
should be reduced, and the images look like 'lips' rather 
than rings. The X-ray images, which appear as rings, indi- 
cate that disordered field may be comparable to the mean 
field. 

Polarization observation of the nebula provides the 
mean field structure and degree of randomness of the field. 
The polarization of the Crab Nebula is observed in ra- 
dio bands (Wilson 1972, Velusamy 1985), in optical bands 
(Oort & Walraven 1956, Schdmidt et al. 1979, Hickson & 
van der Bergh 1990, and Michel et al. 1991), and in X-ray 
bands (Weisskopf et al. 1978). The X-ray polarization is not 
spatially-resolved but a mean polarization is obtained to be 
~ 20%. Almost the same value is obtained also in optical 
bands (Oort & Walraven 1956). 

In this paper, we calculate polarization properties of 
a nebula which expands with relativistic speeds and with 
disordered magnetic field as well as mean toroidal fields. We 
then compare the results with the Crab observations, so that 
the magnitude of the disordered magnetic field is estimated. 



2 POLARIZATION MODEL 

2.1 Polarization from a medium in relativistic 
motion 

Let us consider synchrotron radiation from a relativistic 
plasma moving relative to the observer at the velocity V. 
The volume emissivities in terms of the Stokes parameters 
in the 'flow frame' comoving with the plasma may be given 
by (Rybicki k Lightman 1979) 

f dl'/ds 
dQ'/ds 
dU'/ds 
\ dV'/ds 
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B' — \B'\ is the magnetic field strength, V is the gamma 
function, — e,m are the charge and the rest mass of an elec- 
tron, and c is the speed of light. Here, we have assumed for 
a given volume element in which the magnetic field B' is 
uniform, and the relativistic electrons and positrons have a 
power- law energy distribution, f{^) — K^~^ , where 7 is the 
Lorentz factor of the particles; p and K are functions of posi- 
tion. The pitch angle distribution is assumed to be uniform. 
In the flow frame, the magnetic field B' and the unit vector 
toward the observer make the angle 0' , and the angular 
frequency of the radiation observed at co is Doppler-shifted 
to a;^ 

If the ideal-MHD condition E + P x B = holds, 
where E and B are the electric and magnetic fields in the 
'observers frame', and /3 = V /c, then the external electric 
field in the flow frame disappears. In this case, we can use 
the argument by Bjornson (1982) that 'the observed position 
angles' ^ in the flow frame and the observer frame are the 
same (see Fig. [1]). The Stokes parameters for the observer 
follow (Blandford & Konigl 1979, Bjornson 1982, Lyutikov 
et al. 2003) 
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where T/ = (1 - /3^)"^/^ and V = - /3cosa)"^ are 

the Lorentz factor and the beam factor of the flow, respec- 
tively, a is the angle between the direction of the observer 
n and the flow velocity V, xo is the offset caused by the 
practical position angle which is measured from the north 
on the sky for a given observer, and uj' — uj/V. Numeri- 
cal integration along the line of sight will be done by using 
grid cells constracted in the observer's frame (see Section 
2.2). Quantities in (3) such as P, xo, «f and uo' are different 
in each grid cell because flow velocities are different cell by 
cell. Suitable transformations are applied for each cell. To 
each cell, the orthogonal coordinates are constructed by the 
unit vectors V = V /\V\ in the x-direction, x V in the 
^-direction and V x { n' x V ) in the z-direction, respec- 
tively. The observed position angle ^ is measured from the 
projected z-axis on the sky, i.e.. 



tan g = — — - cos a + — - sm a , 



By 
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where 5^, By, B'^ are the magnetic field components in the 
flow frame. The offset angle xo is defined by the equation 
Xo = — cos""*^ (Z • 7>"), where 7." = — sinaV + cosaz and Z 
is the unit vector of the north direction on the sky. In Pa- 
per I, we used the transformation of the volume emissivity 
with the factor F/P^, and this is not correct. The factor 
should be used. Fig. 2 (bottom) of Paper I may be re- 
placed by Fig. [2ja) of the present paper. However, there is 
no apparent difference between them; the lip shape persists. 
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Figure 1. (a) and (b) show the flow frame and the observer 
frame, respectively. The position angle ^ is an invariant. The unit 
vectors of the flow velocity, the observer's direction and the elec- 
tric vector of the linearly polarized wave are indicated by V, n 
and e, respectively. 



2.2 Construction of polarization map 

If the nebula is optically thin, it is straightforward to ob- 
tain the observed Stokes parameters, simply by integration 
of (|3|) along the line of sight. We use the observer's frame 
(X,y,Z) such that the observer locates at X = +00, the YZ 
plane defines the sky with the north in Z-direction, and the 
nebula center (pulsar) is located at the origin. The volume 
emissivities of (|3|) are functions of X,Y,Z and n, so that 
maps of the stokes parameters are given by 



QAY,Z) = J 
U.{Y,Z) = J 



dIu;{X,Y,Z,n) 



dX 



dX, 



+ CXD 



+ 00 
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We construct grid points on the plane of sky (Yi^Zj) for 
polarization maps, and on the line of sight Xk for integration 
([5]) with increment of As. The Stokes parameters as the 



whole nebula are obtained by 

i 3 
i 3 

Ut = ^^U{Y,,Z,), 



(6) 



with which the mean polarization degree P and position 
angle (P.A.), x, are obtained by P = \/ Qt^ + Ut^ / It and 
X = (l/2)tan~^ (Ut/Qt), respectively. 



2.3 Crab Nebula model 

We start our calculation with a simple disc model with pure 
toroidal magnetic field (Paper I). The disc represents the 
post shock flow with an inner radius Rs located at the shock 
and with a constant semi-opening angle Oq. The flow is based 
on KC: the radial flow velocity V{R) is given (KC), where 
= (X^ + + Z^)^/^. The post shock flow is character- 
ized by the a parameter which is the ratio of Foynting flux 
to the kinetic energy flux just before the shock. The nebula 
flow suffers adiabatic and synchrotron losses. In Paper I, we 
calculate the evolution of distribution function. The param- 
eters K and p in ([T]) are determined as functions of R so as 
to fit the distribution function obtained in Paper I at each 
point in the nebula. The magnetic field distribution Bkc (R) 
is also given by the KC model. We use following a parameter 
set: a — 0.003, the wind luminosity = 5 x 10^^ erg s~^, 
the wind Lorentz factor — 3 x 10^, the shock distance 
= 3 X 10^^ cm, the power-law index at the shock ps — 3, 
and the thickness of the disc — ±10° degree. The inclina- 
tion angle of the axis of the disc to the observer \s i — 28° 
(Weisskopf et al. 2000). 



2.4 The case of disordered toroidal field 

It has been suggested in Paper I that the nebula field is not 
'pure toroidal' but may be dominated by disordered fields. 
Scale length of the randomness (turbulent spectrum) is not 
known; it can be microscopic or just below the resolution 
of observation. In this section, let us consider the case of 
disordered magnetic field. 

In general, the magnetic field can be decomposed to 
the mean field B'q and the random field B'l. The degree of 
randomness may be characterized by 



< -B'l > 
5^2+ < B'l^ > ' 



(7) 



where <> indicates the spatial average in a scale below the 
observational resolution and the primes indicate that h is 
evaluated in the flow frame (< B'l^ > is still a function of 
position). The synchrotron radiation for such cases is stud- 
ied by Korchakov k Syrovat-skii (KS) (1962). From (3), in- 
cluding relativistic motion of the flow, we have the observed 
Stokes parameters: 
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where 



(a) Case KC Intensity 



(b) Case KC Polanization degree 



"As 



s + As 



B\ 1 



(9) 



, (10) 



B' = VIST = VBi)"" + 2Bo • B'l + 



cos 2x = 2 



X n'l \B' X n'l 



(11) 
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The volume emissivities are obtained for each grid cell 
and are integrated numerically according to (O . Given by a 
nebula model, macroscopic quantities such as V and Bo are 
assumed to be constant in a given grid cell, but gradually 
change cell by cell. In each cell, averaging for the random 
field Q and ([TO]) is done by Monte Carlo method. KS use 
the special coordinate in which U vanishes. By using this 
coordinate, only Q remains, and Monte Carlo integration 
can be only one time; the factor cos 2% appears in (p^ . but 
does not sin 2%, where x is defined by ()12p and represents the 
angle between the polarization E- vector and the projected 
magnetic field. The value of x varies place by place within 
the cells. Redistribution of the Stokes parameters for the 
observers coordinate can be done by dS]) by using the angle 
^0 defined by the projected mean field: 



tan ^0 = — ^7— cos a + — — sm a . 



(13) 



J, and that the plasma's distribution functions are 



As stated above, integration along the line of sight has 
been done in use of Monte Carlo method with randomly 
distributed B'lS. In the present calculation, we assume for 
comparison purposes that the energy density of the nebula 
field is equal to the KC value B]^c, i.e 

the same as Paper I. It is obvious that the disordered fields 
change flow dynamics: (i) the disordered fields contribute 
to the magnetic pressure rather than magnetic tension, (ii) 
formation of the disordered fields must be associated by 
heating, and (iii) the distribution function must be changed. 
These issues are postponed to subsequent papers. 

Apart from the disordered field, it is also suggested by 
some authors (Hester et al. 2002, Mori et al. 2004, Ng & 
Romani 2004, and Paper I) that if the front-back (north- 
west to south-east) intensity contrast is due to Doppler 
boost, the flow velocity is much faster than the KC flow 

2000 km s~^) and will be ^ 0.2c. Therefore, we also 
calculate the case with much faster flows. 

In summary, we calculate three cases; Case KC, in which 
plasma properties follows the KC model (see Section 2.3 in 
detail); Case D, in which we include the disordered field 
with b (the degree of randomness) . For comparison purposes 
the total magnetic energy density and flow speed are left 
unchanged; Case DR, in which, in addition to Case D we 
change the flow velocity to be 0.2c through out the whole 
nebula by hand (we ignore the flow dynamics) . 



(c) Case KC P.A.(E-vector) 



Polarization degree (% 



(d) Case KC Polarization degree along the ring 




270 360 



Figure 2. Polarization properties of the pulsar disc nebula 
for Case KC, observed at 5.2 keV. (a)Intensity map. The unit 
length is a shock radius and the unit of the intensity is 0.016 
erg s"-*^ cm~^ str"-*^ eV"-*^. (b) Polarization degree map drawn only 
for the regions where the intensity is strong enough, regarding cal- 
culation errors. (c)P.A. map (E- vector), (d) Polarization degree as 
a function of the azimuthal angle along an ellipse with the semi- 
major axis of 2.5Rs and semi-minor axis of 1.5Rs. The labels 
'SW, 'NE' and 'SW below the horizontal axis indicate south- 
west, north-east and south-west, respectively. 



3 RESULTS 

3.1 Case of Pure Toroidal Field 

Fig.[2]shows the results of Case KC. The observation energy 
is assumed to be 5.2 keV. The intensity map (Fig. EJa)) 
shows a lip-shaped nebula, as was pointed out in Paper I. On 
the major axis (NE-SW) of the ring (perpendicular to the 
rotation axis of the pulsar) and its vicinity, the pitch angles 
of the particles radiating toward the observer are small, and 
as a result synchrotron emissivity is decreased. The dimness 
is cased not only by decrease of single-particle emissivity but 
also by decrease of number of the particles emitting to the 
observing bands. Fig. EJb) shows the polarization degree. In 
Fig. [2jd) , the polarization degree is plotted as an azimuthal 
function along the ring. On the minor axis (rotational axis 
of the pulsar) , the polarization degree reaches the maximum 
value of ^ 75%. There is depolarized regions along the major 
axis, where the polarization degree decreases to ^ 55%. 

In Case KC, the mean polarization degree as the whole 
nebula is 52.3%, and P.A. is found to be 132° at 5.2 keV. 
This polarization degree is much larger than the observed 
value (- 20%) in X-ray (Weisskopf et al. 1978). Although 
we calculate also for optical bands, the results are hardly 
changed so that they are not shown in figures. Since the 
mean polarization in optical is ^ 19% (Oort and Walraven 
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(b) Case D Polarization degree 
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Figure 3. Polarization of the entire nebula as functions of the 
inclination angle i between the axis of the disc and the line of 
sight for b = 0.0, 0.2, 0.4, 0.6 and 0.8. The empirical curves P = 
75.8(1 - 6)(cosz)2-^^ are overlaid. 



1956), the calculated polarization degree in optical is again 
much higher than the observed. 



(c) Case D P.A.(E-vector) 
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(d) Case D Polarization degree along the ring 
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3.2 Effects of Disordered Field and Relativistic 
Motion 



Figure 4. Same as Figure [2] but for Case D with b = 0.6. 



In Case D, we change the parameter b to find the value 
which reproduces the observed polarization degree. Fig. [3] 
shows how the mean polarization behaves as a function of b 
and the inclination angle i, the angle of the symmetry axis 
of the nebula from n (the line of sight). The result can be 
fitted by an empirical relation, P = 75.8x(l — 6)x (cosi)^"^^. 
For the Crab Nebula, the inclination angle is known to be 
i = 28° (Weisskopf et al. 2000), and therefore the observed 
polarization degree can be reproduced if 6 = 0.6. 

Polarization properties for Case D with b — 0.6 are sum- 
marized in Fig. Uat the observation energy of 5.2 keV. As 
has been mentioned in Section 1, almost the same polariza- 
tion degrees are observed both in X-ray and optical bands. 
The calculated polarization properties for the optical bands 
are hardly different from those for the 5.2 keV band. This 
is thus consistent with the observation. The polarization de- 
gree takes the maximum value of ^ 40% on the northwest 
side and the southeast side of the ring. This also agrees with 
the observation in optical bands (Schmidt et al. 1979, Hick- 
son & van der Bergh 1990, Michel et al. 1991). 

Case D is favored by at least two observational facts 
that (1) the intensity map appears in a ring (not lip-shaped), 
and that (2) the polarization degree is ~ 40% at the highest 
and ~ 20% on average. This indicates that the disordered 
magnetic component is 60% of the total. 

The front-back contrast in Case D is smaller than the 
observed value 3.4 (Mori et al. 2004) because the model fol- 
lows the KC flow. In Case DR, the flow velocity is assumed 
to be 0.2c in the whole nebula, and thereby front-back con- 
trast becomes 3.1 (Fig.[5|). The average polarization of Case 
DR is increased only by 0.14% as compared with Case D. 
This increase is caused by the high polarization degree at 
the north-west region weighted by Doppler boosting. Rel- 
ativistic effects is rather prominent in distribution of the 



depolarization. When the nebula flow is non-relativist ic as 
seen in Case KC and Case D, the depolarization takes place 
along the major axis symmetrically, but in Case DR the de- 
polarized regions shift slightly to the northwest side which 
expands toward us (Fig. [5jb) and (d)). The electric field 
vectors of the synchrotron radiation shown in Fig. [5jc) is 
no longer vertical to the projected magnetic field on the 
sky as the flow becomes relativistic. These relativistic effect 
is difficult to see in the complex optical polarization maps. 
However, it might be verified if spatially-resolved polariza- 
tion in X-ray bands is performed because the ring structure 
has a better symmetry in X-ray, and obscuring by filaments 
is less. After all, the relativistic effect does not change the 
estimate of the disordered magnetic field. 

We reproduced the maps of Stokes parameters, Q-map 
and U-map in Fig. [6l which are compared with the obser- 
vation (Michel et al. 1991). Because interaction of the fil- 
aments and the nebula flow changes the direction of the 
magnetic field, comparison may be limited within the area 
of l' around the pulsar. The shape of hourglass (Michel et 
al. 1991) is well reproduced in our Q map. It follows from 
this fact that the average magnetic field is toroidal. Because 
depolarization occurs along the major axis on which U has 
positive value, the absolute value of U is small (see Fig. [6]). 
Therefore, U-map is easy to be influenced by disturbances, 
and therefore comparison of the observations to the model 
is difficult to make. 



4 SUMMARY AND DISCUSSION 

We calculate the synchrotron polarization properties from a 
relativistically expanding disk which has both the toroidal 
mean field component and the disordered component. We 
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Figure 5. Same as Figure [2] but for Case DR with b = 0.6, 
Vr = 0.2c. (d) shows the comparison between the polarization 
degree of Case D (soHd hue) and that of Case DR (dotted hne). 



Figure 6. The comparison between the Q, U maps observed in 
optical bands (Michel et al. 1991) and those of Case DR (b = 
0.6, [3 = 0.2). The upper row is an observation, and the lower is 
the model. The Q map (left) is N/S polarized light minus E/W 
polarized light, the U map (right) is NW/SE polarized light minus 
NE/SW polarized light. 



compared our results with observations of the Crab Neb- 
ula. The distribution function of plasma particles and the 
total magnetic energy density (Bq + Bf)/87T in the disk are 
given as a function of the distance according the KC 
flow with a = 0.003. The best fit model reproducing the ob- 
served polarization indicates that the disordered magnetic 
field contributes 60% of the total magnetic energy. This re- 
sult is also supported by the fact that if the magnetic field 
is pure toroidal, then the nebula image becomes lip-shaped, 
while if there is such a disordered magnetic field, it becomes 
ring-shaped as observed. We also attempt to reproduce the 
front-back contrast of the ring by changing the flow velocity. 
The estimate of disordered field is not changed even if a high 
speed flow of 0.2c is introduced. It is shown that relativist ic 
motion of the expanding nebula can be detected as asym- 
metry in P.A. and in distribution of depolarized regions. 

Nothing is known about wavelengths of the disordered 
magnetic field according to our study. It is pointed out that 
series of magnetic neutral sheets in the wind can be col- 
lapsed at the shock (Lyubarsky 2003). If it is so, the disor- 
dered magnetic field may be formed in the process of mag- 
netic reconnection. The scale lengths of the disordered field 
is thought of as c/^l ^ 1.6 X 10^ or less. 

Recently, relativistic magneto-hydrodynamic (RMHD) 
simulations are done for the flow after the termination shock 
of the pulsar wind, and formation mechanism of the disk 
and the jets is discussed (Komissarov & Lyubarsky 2003, 
Del Zanna et al. 2004). Because of the axisymmetry, the 
magnetic field is assumed to be pure toroidal in those sim- 
ulations. Therefore, even though vortex is formed in the 



meridional plane, it does not cause significant depolariza- 
tion (Bucciantini et al. 2005, Del Zanna et al. 2006). 

It is notable that if the disordered field dominates, then 
the magnetic field contributes more as pressure than ten- 
sion in the flow dynamics, as compared with the case of 
pure toroidal field, which is assumed in the previous RMHD 
simulations. As a result, magnetic pinch effect may be over- 
estimated in the previous simulations. If the disordered mag- 
netic field is taken into account, the jets may be weaken, and 
magnetic pressure may cause acceleration of the disk flow. 

Hickson (1990) finds highest polarization of ^ 60% in 
a region 1^8 apart from the pulsar in a south-south-east 
direction and the eastern bay. For these regions, very ordered 
magnetic field is expected. However, polarization degrees of 
35% for the wispl, 23% for the wisp3 (Hickson 1990), and 
~ 30% of the wisp3 (Scargle 1971) seem to be consistent 
with b ^ 0.6. 

For a future possibility, we simulated change of polariza- 
tion for the entire Crab Nebula through a lunar occultation. 
Fig. [71 shows the results when the nebula is gradually hid- 
den, while Fig.[8]shows a recovery phase. The direction of the 
occultation can be seen in the images of the figures. In the 
disappearance phase (Fig. [7]), the total polarization degree 
decreases at first because of hiding highly polarized parts, 
and then it is minimized and increases. The late phase will 
not be observed because the south-east part of the ring is 
dim, and the radiation from the jet will dominate. In the 
reappearance phase, the polarization degree start with a 
very high state ~ 40% and gradually recovers. There will 
be no significant change of position angle. 
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(b) (c) 

Figure 7. Polarization at a lunar occultation of the Crab Neb- 
ula (disappearance phase). We use Case DR, b = 0.6, f3 = 0.2. 
(a)Appearance of luna occultation from the northwest to the 
southeast along symmetrical axis of the nebula, (b) Change of 
averaged polarization degree. (c)Change of P.A. The labels 0~7 
of the horizontal axis of panel (b) and (c) correspond to eight 
figures from on the upper left to lower right of sub-panels in (a). 




(b) (c) 

Figure 8. The same as Fig. [71 but for the appearance phase. The 
labels 1~8 of the horizontal axis of panel (b) and (c) correspond 
to eight figures from on the upper left to lower right of sub-panels 
in (a). 
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